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^ . Abstract 
OXJ, 

This paper shows a novel calculation of the mean square displacement of a classical Brownian 
particle in a relativistic thermal bath. The result is compared with the expressions obtained 
by other authors. Also, the thermodynamic properties of a non-degenerate simple relativistic 
gas are reviewed in terms of a treatment performed in velocity space. 



1 Introduction 
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, Interest in relativistic statistical mechanics has been revived due to its applicability in several fields 

' of study and the need of its internal theoretical consistence |lj[2]. New experimental devices, such 

, as relativistic heavy ion coUiders, have also drawn attention to this topic since they provide a 

new framework to test different approaches to relativistic transport theory [3]. In the context of 
stochastic processes theory, formalisms involving relativistic thermodynamics have been proposed in 
' order to understand the nature of how relativity may broaden the usual conceptions of mathematical 

probability theory. 

In this sense, the analysis of processes such as Brownian motion in a relativistic framework 
deserves a deeper look. Other issues related to the inclusion of non-equilibrium effects in relativity 
theory have been addressed in previous work [4] [5]. In the case of large (relativistic) temperatures, 
a description involving a Juttner function is required in order to compute the mean square dis- 
placement of Brownian particles [6]. This calculation is described in the present work. The interest 
of Brownian motion has been recently revived due to its applicability in different areas, including 
astrophysics [9j. 

This paper is divided as follows. In section two, the basic properties of a non-degenerate 
relativistic gas are reviewed and compared with their Maxwell-Boltzmann function counterparts at 
mildly and ultra- relativistic regimes. Section three is devoted to the study of relativistic Brownian 
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motion modeling the dynamics of a Newtonian brownian particle in a relativistic thermal bath. Due 
to this aproximation, only elementary mathematical tools are involved in the calculations. Finally, 
some theoretical remarks on the use of relativistic statistical mechanics are included in section four. 



2 Distribution functions 

Some thermodynamic systems cannot be modeled using the usual Maxwell-Boltzmann distribution. 
Instead, a more appropriate equilibrium distribution function should be considered in high temper- 
ature cases. For example, the Jiittner distribution function for a special relativistic non-degenerate 
simple gas at rest is [1]: 

_ 21 

where 

-1 

is the usual Lorentz factor for molecular speed uj, z — is the relativistic parameter at temper- 
ature T for the system consisting of particles of mass to, K2 is the modified Bessel function of the 
second kind, fc and c are Boltzmann's constant and the speed of light, respectively. 

Figure 1 shows a comparison between both speed distribution functions as functions of /3 = ^ , 
corresponding to Eq. ([l]) and a Maxwell-Boltzmann distribution written in terms of 7: 

for an electron gas. Although the differences seem small, they are measurable and can be detected 
in certain astrophysical processes such as the Sunyaev-Zel'dovich effect [7]. 
In special relativity, the four-vector v'' is related to the three velocity w by: 



(4) 
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and the volume element in velocity space reads d^w = 47rw^7^dw [8]. One fundamental difference 
between both distributions, Maxwellian and Jiittnerian, is that the former allows for molecular 
speeds greater than c, while the latter has a cutoff at /3 = 1. Both behaviors are shown in Fig. 1. 
On the first plot both distributions have a similar qualitative behavior, where most particles lie on 
the right side of the peaks. However, for larger values of z, as can be seen in the second plot, both 
distributions are quite different. Notice how the dotted, classical, distribution exceeds the physical 
Hmit = 1 while the Jiittner distribution function vanishes exactly at that value. 

We are now in position to address a simple application of relativistic thermodynamics and the 
corresponding comparison with the non-relativistic case. 



3 A simple approach to relativistic Brownian motion 

Before addressing the calculation of the mean internal energy per particle for a non-degenerate 
gas, both in the non-relativistic and relativistic cases. In the non-relativistic case, this quantity is 
identified with the average integral: 
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Figure 1: This figure shows the probabiHty distribution functions, for z = 0.01 (left) and z = 0.1 
(right), both in the classical case (dotted line) and in the relativistic case (solid line). The small 
differences at z = 0.01, corresponding to an electron temperature of 6 x 10^ K for an electron gas, 
are large enough to account for relativistic corrections to the Sunyaev-Zel'dovich effect in clusters of 
galaxies. The relativistic features of the probability function are significantly enhanced for z = 0.1, 
corresponding to an electron temperature of 6 x lO^K. 
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<e>=y -mv''fMBiv)dv = -m^ \^—] / v^e dv. (5) 
The usual change of variable u = \/-^BpV leads to the result: 
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This result is identified with the energy equipartition theorem. Each degree of freedom contributes 
units of energy to the total internal energy per particle of the system. As we shall see below, 

the energy equipartition theorem does not hold in relativistic thermodynamics. 

In the relativistic case, since the kinetic energy per particle is mc^j, the counterpart of Eq.|[5|) 

reads: 



< e >= / mc^-/f{w){ATrw'^j^)dw. (7) 

In terms of 7 Eq. ((?]) reads: 



zK2 {-) Jl 

which yields an expression for < e > in terms of modified Bessel functions of the second kind, 
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namely [Ij 

<£>^3fcr + mc^ (9 

Equation ([9]) reflects the fact that, for low enough temperatures, internal energy is well described 
by a linear function of the temperature T, and a non-linear behavior is exhibited for increasing z. 
The use of asymptotic expansions for the modified Bessel functions present in this equation, for 
z < 1 leads to the result of Eq. © [lO]. 

In order to account for the mean square displacement < > of a classical Brownian particle of 
mass M in a relativistic gas, we assume that the relativistic parameter z of the gas is not negligible, 
while the equation of motion describing the displacement still reads: 

dv -* 

M—+llV^Fa, (10) 

where Fa is a random force, /i is the dynamic friction coefficient and ~v' is the velocity of the particle. 
As usual, we assume that < Fa >= so that 

f|^')H.). (11) 

In Eq. (fTTj) we will now use the expression given by ^ , for < e > , since the relativistic parameter 
is significant. We finally get, following the usual steps to compute the mean-square displacement 
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t, (12) 



which reduces to the well-known Einstein-Smoluchowsky equation for the non-relativistic case, if 
we substract the rest energy to the internal energy given by Eq. ^ . Historically, this equation has 
been used to establish the value of Boltzmann's constant k from direct measurement of the mean 
square displacement in various experimental situations [9]. In figure 2, we plot the next equation 
in order to compare the relativistic mean square displacement with its non-relativistic counterpar: 

r2 1 ifi(i) 1 

(— )=2(1 + --^--), (13) 
2Dt oz K2{^) oz 

where D = For z <C 1 we obtain the correct non-relativistic limit 



4 Final remarks 

The framework here presented provides a direct prediction for the value of the mean square dis- 
placement of a massive particle exhibiting Brownian motion in a mildly relativistic, simple gas. 
There is an elementary physical argument to support the approach here proposed. Brownian mo- 
tion is caused by random coUisions of fight particles with massive particles. In a certain range 
of temperatures, the relativistic parameter Z = for the massive particles can be considered 
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Figure 2: Comparison of the values of the mean square displacement for relativistic and non- 
relativistic temperatures. For different values of z the asymptotic behavior of the ratio tends to 
1. 
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negligible, allowing a classical description for their dynamics, while for the thermal bath, z ~ 1 
is significant. Since a non-relativistic distribution predicts an important fraction of particles with 
supraluminal velocities for large z, a relativistic treatment describing the properties of the gas is 
certainly needed. Other calculations of the mean square displacement in this physical situation 
may require the use of numerical methods or more elaborated approaches [T2j . 

A similar calculation has been permormed in reference [12j in the case of a Wiener process in 
which is time-dependent for a given value of z. Work in this direction will be permormed in the 
future. 
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